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THEORETICAL ASPECTS OF RARE KAON DECAYS a 



David Greynat and Eduardo de Rafael 
CPT, CNRS-Luminy, Marseille, France 



Most of the analytic approaches which are used at present to understand Kaon decays, get their 
inspiration from QCD in the limit of a large number of colours N c . After a general overview, 
we illustrate this with a recent application to the evaluation of the process Kl — * i n 
the Standard Model. 



1 INTRODUCTION 

In the Standard Model, the electroweak interactions of hadrons at very low energies are con- 
veniently described by an effective chiral Lagrangian, which has as active degrees of freedom 
the low-lying SU(3) octet of pseudoscalar particles, plus leptons and photons. The underlying 
theory is a SU(3)c x SU(2)l x U{\)w g au g e theory which is formulated in terms of quarks, 
gluons and leptons, together with the massive gauge fields of the electroweak interactions and 
the hitherto unobserved Higgs particle. Going from the underlying Lagrangian to the effective 
chiral Lagrangian is a typical renormalization group problem. It has been possible to integrate 
the heavy degrees of freedom of the underlying theory, in the presence of the strong interactions, 
perturbatively, thanks to the asymptotic freedom property of the SU (3)-QCD sector of the the- 
ory. This brings us down to an effective field theory which consists of the QCD Lagrangian with 
the u, d, s quarks still active, plus a string of four quark operators and mixed quark-lepton oper- 
ators, modulated by coefficients which are functions of the masses of the fields which have been 
integrated out and the scale \i of whatever renormalization scheme has been used to carry out 
this integration. We are still left with the evolution from this effective field theory, appropriate 
at intermediate scales of the order of a few GeV, down to an effective Lagrangian description 
in terms of the low-lying pseudoscalar particles which are the Goldstone modes associated with 
the spontaneous symmetry breaking of chiral-5C/(3) in the light quark sector. The dynamical 
description of this evolution involves genuinely non-perturbative phenomena. There has been 
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recent progress in approaching this last step, using analytic methods formulated within the con- 
text of QCD in the limit of a large number of colours N c (QGD^), (see e.g. re^ for a recent 
review.) 

The strong and electroweak interactions of the Goldstone modes at very low energies are 
described by an effective Lagrangian which has terms with an increasing number of derivatives 
(and quark masses if explicit chiral symmetry breaking is taken into account.) Typical terms of 
the chiral Lagrangian are 

£ eff = J F o tr[D fl UD^)+L 10 tx(u^F Rflu UF^)+--- (1) 
> , ' , ' 

7T7T — *7T7T , K — >7TeV ^ 

+ e 2 Cti(Q R UQ L U^) +••• (2) 
v * ' 

-e 2 C*-^ (n+n~+K+K- ) 

- ^V ud V: s gsF* (D,UD^) ^ + • • • , (3) 
s , ' 

K — >7T7T , K — >7r7T7r 

where U is a 3 x 3 unitary matrix in flavour space which collects the Goldstone fields and which 
under chiral rotations transforms as U — > VrUV^; D^U denotes the covariant derivative in the 
presence of external vector and axial- vector sources. The first term in the first line is the lowest 
order term in the sector of the strong interactions^, Fq is the pion-decay coupling constant in 
the chiral limit where the light quark m asse s u, d, s are neglected (Fq ~ 90 MeV); the second 
term shows one of the couplings at C(p 4 )^El the second line shows the lowest order term which 
appears when photons and Z's are integrated out (Ql = Qr = Q = diag.[(2/3, —1/3, —1/3]), in 
the presence of the strong interactions ; the third line shows one of the lowest order terms in the 
sector of the weak interactions. The typical physical processes to which each term contributes are 
indicated under the braces. Each term is modulated by a coupling constant: Fq, L\q,... C...g%..., 
which encodes the underlying dynamics responsible for the appearance of the corresponding 
effective term. The evaluation of these couplings from the underlying theory is the question we 
are interested in. The coupling g% for example, governs the strength of the dominant AI =1/2 
transitions for K decays to leading order in the chiral expansion. 

There are two crucial observations to be made concerning the relation of these low energy 
constants to the underlying theory. The low-energy constants of the Strong Lagrangian, like 
Fq and L 10 , are the coefficients of the Taylor expansion of appropriate QCD Green's Functions. 
By contrast, the low-energy constants of the Electroweak Lagrangian, like e.g. C and g$, are 
integrals of appropriate QCD Green's Functions. Their evaluation appears to be, a priori, 
quite a formidable task because they require the knowledge of Green's functions at all values 
of the euclidean momenta; i.e. they require a precise matching of the short-distance and the 
long-distance contributions to the underlying Green's functions. These two observations are 
generic in the case of the Standard Model, independently of the 1 /N c expansion. The large-iV c 
approximation helps, however, because it restricts the analytic structure of the Green's functions 
to be meromorphic functions; i.e., they only have poles as singularities. 

The QCDqo approach that we are proposing in order to compute a specific coupling of the 
chiral electroweak Lagrangian consists of the following steps: 

1. Identification of the relevant Green's functions. 



2. Evaluation of the short-distance behaviour and the long-distance behaviour of the relevant 
Green's functions. 



3. Hadronic approximation of the underlying Green's functions with a finite number of poles; 
i.e., the minimum number required to satisfy the leading power fall-off at short-distances, 
as well as the appropriate xPT long-distance constraints. 

We have checked this approach with the calculation of a few low-energy observables: 

i) The electroweak Am n mass difference^. 

ii) The hadronic vacuum polarization contribution to the anomalous magnetic moment of 



the muon 
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Hi) The 7T° — * e + e~ and n — * /i + fi~ decay rates^. 



These successful tests have encouraged us to pursue a systematic analysis of -fT-physics 
observables within the same large-iV c framework. So far, the following calculations have 
been made: 

• The Bk Factor in the Chiral Limi^. 

• Weak Matrix Elements of the Electroweak Penguin Operators Qj a ndCh^ ^. 

• Hadronic Light-by-Light Scattering Contribution to the Muon g — 2^212] 

• Electroweak Hadronic Contributions to the Muon g — 2^3. 

• Weak Matrix Elements of the Strong Penguin Operators Q4 and qJ^. 

It appears now to be possible to apply the same techniques to the study of rare K 
decays as we ll, a nd we shall illustrate this here, with a first application to the process 

k l -+ vL+n-m 



2 RARE KAON DECAYS IN X PT 

In fact, it was at the Blois Conference in '98 where some of the first applications of xPT to 
rare Kaon decays were reported. The predicted invariant mass spectrum of the two gammas 
in the mode Kl — > vr°77, to lowest order in the chiral expansion LLULLiil was discussed prior 
to the earlier experimental measurements by the NA31 and E731 collaborations. Since then, 
%PT has been confirmed as the appr opr iate framework to study rare Kaon decays with many 
successful applications, (see e.g., ref.*^ for the latest conference review on the subject, where 
further references can be found.) 

Unfortunately, the predictive power of xPT is seriously limited at present, because of the fact 
that several coupling constants of the higher order terms in the electroweak chiral Lagrangian 
remain unknown. As a result, many predictions at present rely on models which are not clearly 
related to the underlying Standard Model theory. The analytic approach which we have briefly 
summarized in the introduction offers, however, an interesting possibility to make progress in 
this field. 

Roughly speaking, from a theoretical point of view, there are three types of if-decay modes: 

• The golden modes of xPT 

They correspond to processes which are fixed by xPT only. The chiral loops are finite and 
there are no unknown counterterms at the order one is working in the chiral expansion. 
Examples of that are K$ — ► 77 and Kl — > 7r°77 

• The short-distance golden modes 

They directly test the structure of the Wilson coefficients, because the relevant effective 
weak Hamiltonian for these processes is of the unique type: quark-current x lepton cur- 
rent. The Wilson coefficients have been calculated perturbatively with a sufficient degree 
of accuracy which makes some of these modes an excellent laboratory of possible Physics 



beyond the Standard Model. The paradigm mode of this type is Kl — ► ir®vv (see e.g. 
ref.H^ for a recent review on this mode) 

• The mixed modes 

This is, unfortunately, the largest class. Some modes are sensitive to short-distance 
Physics, CP violation in particular, like Kl — ► 7r°e + e~ and e/e', but their calculation 
in the Standard Model brings in unknown couplings of the chiral Lagrangian. It is here 
where hard theoretical effort is mostly required. 

3 K L -> n+yT 

This is a decay, with both an interesting short-distance component which is particularly sensitive 
to Vtdi and a long-distance contribution which has a large absorptive component from the 
dominant 77 intermediate state 1 ^, illustrated in Fig. 1. The present experimental rates are 

Br(K L -> = (7.18 ± 0.17) x 1(T 9 BNL-E871ED (4) 

Bt(K l -> e+e") = (8.7 J x 10~ 12 BNL-E871^ (5) 




From a phenomenological point of view, it is convenient to normalize the Kl — » // decay rate to 
the one of the Kl — ► 77, which is also known experimentally [Br (if^ — » 77) = (5.86 ± 0.17) x 
10" 4 ]. Then 
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with 

In fact, the contribution from the absorptive amplitude Im^^ almost saturates the observed 
experimental rate, leaving very little room for the real part 

\ReA^\ 2 < (7.1 ± 0.2) x 10~ 9 . (8) 



The real part of An has a component from the 77 loop which is divergent and requires a local 
long-distance counterterm. The relevant couplings from the long-distance effective Hamiltonian 
are the following 

^eff' = C^j-e, va ^A^K L ) + ffl ^>{d,K L ) 1^ . (9) 

The first term provides the if£77~vertex interaction in the loop diagram of Fig. 1, with an 
unknown coupling C which can however be fixed (in modulus only) from the observed Kl — > 77 
decay rate: 

The second term is the long-distance contact term, which appears at the same order in the 
chiral expansion as the contribution from the 77 loop. The coupling Xld(/^) is the quantity that 
one would like to evaluate in the Standard Model. In terms of this coupling, one finds 



Re^ LD >(A) = xM + ^ 
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_ + _ ln( _x,+ c(/3 , ; 



(ii) 



where the function CiBA corresponds to a finite three-point loop integral which is known, and 
can be found in refs. The divergence associated with this diagram has been renormalized 
within the MS minimal subtraction scheme of dimensional regularization. The logarithmic 
dependence on the renormalization scale /x displayed in the above expression is compensated by 
the scale dependence of the renormalized coupling Xld(/-0- Let us stress here that, in contrast 
with the usual situation in the purely mesonic sector, this scale dependence is not suppressed 
in the large-iVc limit, since it does not arise from meson loops. 

The amplitude Aji also gets a tree level contribution from the short-distance Hamiltonian^^ 

wg } = • • • - %--4r ^ Ynl + XtY ^ f^^rM + h.c. , (12) 

en V2 vr sm z &w \ * J 2 

where 
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K = V c * s V cd , \ t = V*V td , x t = ^ r . (13) 

lvl w 

In the Standard Model, this is the term in the effective Hamiltonian which emerges after integrat- 
ing out the heavy degrees of freedom (Z, W; t, b, and c) in the presence of the strong interactions, 
and which couples directly the quark current sj^ d to the lepton current Ij^ ~^ 5 l- Here Y/vl 
and Y(xt) are functions of the masses of the integrated particles, which can be found in refill 
The effect of the short-distance Hamiltonian ()12j) in the amplitude Aji is to induce a shift 
in the Xld (/x)-coupling 

XldO«) =>• Xld(m) - Xsd , (14) 

with x SD fixed by the relation 

- gg- . * Re [X C Y NL + X t Y(x t )} = -Cx XsD . (15) 

Without loss of generality, we can fix \i at the p mass and define the parameter 

Xoh = Xlo(M p ) - Xsd (16) 

The branching ratio TZj, in Eq. © is then a function of only the unknown parameter x c ff> and 
the predicted shapes for I = fj,,e are as in Fig. 2. From the comparison with the experimental 
value of 7Zp,[i (the horizontal band in Fig. 3a) we conclude that 

3.9 < Xoff < 6.5 . (17) 



The minimum of the predicted parabola in Fig. 2a corresponds to the lower bound when 
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Fig. 2a The predicted branching ratio TZp,^ as a function ofXed- The horizontal lines corre- 
spond to the present experimental result with errors. 
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Fig. 2b The predicted branching ratio TZee as a function of Xcb- The horizontal lines corre- 
spond to the present experimental result with errors. 

2 

Because of the small electron mass, the term log j-A dominates the r.h.s. of Eq. and there is 
very little we can learn on x c s from the, otherwise, quite remarkable experimental determination 

Of TZee- 



3.1 EVALUATION OF Xld(^) IN THE LARGE-N C QCD FRAMEWORK 

The relevant Green's function for the process we are concerned with is the four-point function 



W^(q,p) = lim| d'xe^J tfyj™ J d 4 ze il - z (0\T {j^x)J u (y)V(0)C^ s=1 (z)}\0) (19) 

= jjw ( l> \ q+p y_ M 2 K ^[<l^P^ (<1 + P) 2 } , ( 2 0) 

where, in the first line, J M is the electromagnetic current; V = qi^^-q the interpolating 
density current operator of the -ftT^-state and C^f =1 the strangeness changing chiral effective 
Lagrangian. The function W[q 2 ,p 2 , (q + p) 2 ] governs both the Kl — ► (i + fi~ and the Kl — » 77 
decay rates. In particular, the relation to the constant C introduced in Eqs. @ and ()1U|) is as 
follows 

W[0,0,M|] = ^!<||>1 C , (21) 

while the evaluation of the Kl — * amplitude requires the knowledge of the function 

W[q 2 , q 2 , M 2 K \ at all values < -q 2 < 00. 

It is well known that the lowest order evaluation of the constant C in SU(3)l X SU(3)r xPT 
gives a contribution which is proportional to the lowest order Gell-Mann Okubo mass relation: 
3M 2 — AM^ + M 2 = 0. The first non-trivial contribution to C comes from chiral loops and a 
combination of undetermined couplings of the effective C^ s=1 Lagrangian at 0(p A ). This is why 
it has become so difficult to evaluate Kl — ► 77, and hence Kl — ► in the Standard Model. 

We do not attempt to do a calculation of C, but we shall use a crucial observation which, 
to leading non-trivial order in the chiral expansion, relates the residue of the short-distance 
behaviour (— q 2 — ► 00) of the function W[q 2 ,q 2 ,M^] to the constant C, known in modulus from 
the observed Kl — > 77 decay rate b . The relation is the following c 

Mm Win 2 a 2 M 2 ] - ^ (c 4- 8 AT M l M 2 K - M 2 \ 



'/ . 



and is based on the following facts: 



The first non-trivial contribution to the constant C in xPT and to leading order in the 
\/N c expansion, comes from the bosonization of the T-product in Eq. (|19|) . with the 
contribution to J ll {x)J v {y) induced by the local one-Goldstone contribution from the 0(p 4 ) 
Wess-Zumino Lagrangian, the 0(p 2 ) bosonization of V(0) and the problematic 0(p 4 ) 
contribution from £.^f =1 (z) d 



• When performing the operator product expansion of the two electromagnetic currents, 
the leading contribution is governed by the local operator: qQ 2 j^q, which also has to be 
bosonized in terms of one-Goldstone states and produces the same combination of fields 
as the one emerging from the Wess-Zumino J^(x)J u (y) product discussed in the previous 
item. The resulting T-product of this gQ 2 7^ (/-operator, with the remaining V(0) and 
C^f =1 (z) operators, has to be evaluated at long-distances. To first non-trivial order, 
there are two pieces which contribute to its bosonization: the one from qQ 2 ^q at 0{p 2 ) 

b It is in this sense that we differ from the large- N c inspired U(3)l X U(3)r phenomenological analysis of refill 
(See also the discussion irJ^l) 

c N ot ic e the the short-distance behaviour obtained here is in disagreement with the hadronic ansatz proposed 
in ref.EEI. 

d We have identified the combination of 0(p 4 ) couplings which contribute 



with V(0) at 0(p 2 ) and £^f =1 (z) at 0(p 4 ) produces a term proportional to the same 
contribution discussed in the previous item; the other one from qQ 2 j l3 q at 0(p A ), which 
is proportional to the L5 coupling, times V(0) at 0(p 2 ) and C^ s=1 (z) also at 0(p 2 ). 
Fortunately, this second contribution brings in known couplings of the chiral Lagrangian. 

• The other possible contributions to W[0, 0, M^] and lim_ (J 2^ 00 W[q 2 , q 2 , M^], at the same 
order in the chiral expansion, either give vanishing contributions, or are proportional to 
the Gell-Mann Okubo combination of masses. However, higher order contributions in the 
chiral expansion may very likely destroy this dynamical symmetry observed at the first 
non-trivial order. This is why our main result in Eq. (|22|) is only valid to first non-trivial 
order in the chiral expansion. Furthermore, the effect of chiral loops, subleading in the 
1/N C expansion are also neglected. 

The procedure to evaluate the effective coupling x L d(aO is now entirely similar to the one 
discussed in ref.^ with the result 
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1.19 ±0.16 



(23) 



where the first db sign comes from the fact that the constant C is only known in modulus. This 
implies a two-fold result for the effective coupling Xca m Eq. (|16|) which, using the short-distance 
evaluation Xsd = ±(1.8 ± 0.1), results in 



2.33 ±0.9 -(1.8 ±0.1) =0.5 ±0.9 1 
Xeff ~ 2.03 ± 0.9 ± (1.8 ± 0.1) = 3.8 ± 0.9 J 

The corresponding branching ratios for the two solutions are shown in Figure 3 below 



(24) 
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Fig. 3 Branching ratio versus Xeff- The horizontal band corresponds to the experimental 
value with errors . The predicted solutions are the two solid vertical bands with vertical error 
bars. 



The higher solution is in perfect agreement with the experimental determination, the lower 
one is in slight difficulty. 



It is very likely that considerations similar to the ones illustrated here in the case of Kl — ► 
fi + [J.~, may also apply to other rare .fT-decay processes and open, therefore, a new theoretical 
perspective in the field of rare K decays. 
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